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ABSTRACT 

The history of the transition from a neutral intergalactic medium (IGM) to one 
that is almost fully ionized can reveal the character of cosmological ionizing sources. 
We study the evolution of the volume filling factors of H ii and He iii regions in a 
clumpy IGM, and discuss the implications for rival reionization scenarios of the rapid 
decline in the space density of radio-loud quasars and of the large population of star- 
forming galaxies recently observed at z > 3. The hydrogen component in a highly 
inhomogeneous universe is completely reionized when the number of photons emitted 
above 1 ryd in one recombination time equals the mean number of hydrogen atoms. 
If stellar sources are responsible for keeping the IGM ionized at z = 5, the rate of 
star formation at this epoch must be comparable or greater than the one inferred from 
optical observations of galaxies at z ~ 3, and the mean metallicity per baryon in the 
universe ~ 1/500 solar. An early generation of stars in dark matter halos with circular 
velocities Vdrc ^ 50 kms^^, possibly one of the main source of UV photons at high- 
z, could be detectable with the Next Generation Space Telescope. Models in which 
the quasar emissivity declines rapidly at z ^ 3 predict a late He ii reionization epoch, 
a feature that could explain the recent detection of patchy He ii Lya at z = 2.9 by 
Reimers et al. (1997) and the abrupt change observed by Songaila (1998) at about the 
same epoch of the Si iv /C iv ratio, but appear unable to provide the required number 
of hydrogen- ionizing photons at z ~ 5. 



Subject headings: cosmology: theory - diffuse radiation - intergalactic medium - galax- 
ies: evolution - quasars: general 
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1. Introduction 

The existence of a filamentary, low-density intergalactic medium (IGM), which contains the 
bulk of the hydrogen and helium in the universe, is predicted as a product of primordial nucleosyn- 
thesis (e.g., Copi, Schramm, & Turner 1995) and of hierarchical models of gravitational instability 
with "cold dark matter" (CDM) (Cen ct al. 1994; Zhang, Anninos, & Norman 1995; Hernquist 
et al. 1996). The application of the Gunn-Pctcrson (1965) constraint on the amount of smoothly 
distributed neutral material along the line of sight to distant objects requires the hydrogen compo- 
nent of the diffuse IGM to have been highly ionized by z ~ 5 (Schneider, Schmidt, & Gunn 1991), 
and the helium component by z w 2.5 (Davidsen, Kriss, k. Zheng 1996). The plethora of discrete 
absorption systems which give origin to the Lya forest in the spectra of background quasars are 
also inferred to be strongly photoionized. It thus appears that substantial sources of ultraviolet 
photons were present at z > 5, perhaps low-luminosity quasars (Haiman &: Loeb 1998) or a first 
generation of stars in virialized dark matter halos with Tvir ^ 10^ K (Couchman & Rces 1986; 
Tegmark, Silk, & Blanchard 1994; Fukugita & Kawasaki 1994; Ostriker & Gnedin 1996; Ciardi k 
Ferrara 1997; Haiman k Loeb 1997; Miralda-Escude k Rees 1998). Early star formation provides a 
possible explanation for the widespread existence of heavy elements in the IGM (Cowie et al. 1995), 
while rcionization by QSOs may produce a detectable signal in the radio extragalactic background 
at meter wavelengths (Madau, Meiksin, k Rees 1997). Establishing the character of cosmological 
ionizing sources is an efficient way to constrain competing models for structure formation in the 
universe, and to study the collapse and cooling of small mass objects at early epochs. While the 
nature, spectrum, and intensity of the background UV flux which is responsible for maintaining 
the intergalactic gas and the Lya clouds in a highly ionized state at 2 < 3 has been the subject 
of much debate in the last decade (e.g., Bechtold et al. 1987; Shapiro k Giroux 1987; Songaila, 
Cowie, k Lilly 1990; Miralda-Escude k Ostriker 1990; Madau 1992; Sciama 1995; Haardt k Madau 
1996; Giroux k Shapiro 1996), it is only in the past few years that new observations have provided 
reliable information on the presence and physical properties of sources and sinks (due to continuum 
opacities) of UV radiation in the interval 3 < z < 5. 

This is the third paper in a series aimed at a detailed study of the production, absorption, and 
reprocessing of ionizing photons in a clumpy universe. In Paper I (Madau 1995) we showed how the 
stochastic attenuation produced by neutral hydrogen along the line of sight can seriously distort 
our view of objects at cosmological distances. In Paper II (Haardt k Madau 1996) a numerical code 
was developed which follows the propagation of photoionizing radiation throTigh an IGM in which 
the H I and He II region networks have fully percolated at the redshifts of interest. It is the purpose 
of this paper to focus instead on the candidate sources of photoionization at early times and on 
the time-dependent reionization problem, i.e. on the history of the transition from a neutral IGM 
to one that is almost fully ionized. The primary motivation of this study can be found perhaps in 
the simple realization that the breakthrough epoch (when all radiation sources can see each other in 
the Lyman continuum) occurs much later in the universe than the overlap epoch (when individual 
ionized zones become simply connected and every point in space is exposed to ionizing radiation). 
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and that at high redshifts the ionization equihbrium is actually determined by the instantaneous 
UV production rate. The fact that the IGA-I is still optically thick at overlapping, coupled to recent 
observations of a rapid decline in the space density of radio-loud quasars and of a large population of 
star-forming galaxies at z ^ 3, will have some interesting implications for rival ionization scenarios. 

Throughout this paper we will adopt an Einstein-de Sitter universe (go = 0.5) with Hq = 
50/i5o kms-^Mpc-i. 



2. The Universe After Complete Overlapping 

The complete reionization of the universe manifests itself in the absence of a Gunn-Peterson 
absorption trough in the spectra of galaxies and quasars at z < 5. In the presence of a uniform 
medium of H i density nmiz) along the path to a distant object, the optical depth associated with 
resonant Lya scattering at Xa{l + z) is 

^GP(^) = = 8.3x10 (1) 

where H is the Hubble constant, /« is the oscillator strength of the transition, and all other symbols 
have their usual meaning. The same expression for the opacity is also valid in the case of optically 
thin (to Lya scattering) discrete clouds as long as nni is replaced with the average neutral density 
of individual clouds times their volume filling factor. The strongest limit on the amount of diffuse 
intergalactic neutral hydrogen at z = 4.3 is provided by a direct measurement of the quasar Q1202- 
0725 flux in regions of the spectrum where lines are absent (Giallongo et al. 1994): a best estimate 
of tgp = 0.02 ± 0.03 is found. Even if 99% of all the cosmic baryons fragment at these epochs 
into structures that can be identified with QSO absorption systems, with only 1% remaining in a 
smoothly distributed component (cf. Zhang et al. 1998), the implication is a diffuse IGM which is 
ionized to better than 1 part in 10^. Except at the highest column densities, discrete absorbers arc 
also inferred to be strongly photoionized. From QSO absorption studies we also know that neutral 
hydrogen accounts for only a small fraction, ~ 10%, of the nucleosynthetic baryons at early epochs 
(e.g. Lanzetta, Wolfe, & Turnshek 1995). 

In this section we shall discuss the last stage of the reionization process, when individual 
ionized zones have overlapped, the reionization of hydrogen (and helium) in the universe has been 
completed, and the IGM has been exposed everywhere to Lyman continuum photons. An often 
overlooked point is worth remarking here, namely the fact that even if the bulk of the baxyons in 
the universe are fairly well ionized at z < 5, the residual neutral hydrogen still present in the Lya 
forest clouds and the Lyman-limit systems (LLS) along the line of sight significantly attenuates the 
ionizing flux from cosmologically distant sources. In other words, while the complete overlapping of 
I-fronts occurs at z > 5, the universe does not become optically thin to Lyman continuum photons 
until much later, at a "breakthrough" redshift Zbr ~ 1-6 (see § 2.3). It is only after breakthrough 
that every point in space at any redshift z < Zbr will be exposed to ionizing photons from all 



-4- 



radiation sources out to Zbr- 



2.1. Cosmological Radiative Transfer 

The radiative transfer equation in cosmology describes the time evolution of the specific in- 
tensity J of a diffuse radiation field (e.g., Peebles 1993): 

d a d \ ^ a ^ ^ c 
at a Of J a Att 

where a is the scale parameter, c the speed of the light, k the continuum absorption coefficient per 
unit length along the line of sight, and e is the proper space-averaged volume emissivity. The mean 
specific intensity of the radiation background at the observed frequency Vq, as seen by an observer 
at redshift Zo, is then 

where v = fo(l + 2;)/(l + 2;o)) and dl/dz is the line element in a Friedmann cosmology. The effective 
optical depth t^s due to discrete absorption systems is defined, for Poisson-distributed clouds, as 

fZ POO «2 J\T 

.=.(.o..„,.)=£<!.'X ■'^H.g^d-O (4) 

(Paresce, McKee, & Bowyer 1980), where d^N/dNmdz' is the redshift and column density distri- 
bution of absorbers along the line of sight, and r is the Lyman continuum optical depth through 
an individual cloud. 



2.2. Intervening Absorption 

The actual distribution of intervening clouds is still quite uncertain, especially in the range 
-|^q16.3 < < iQi'^-^cm"^, where most of the contribution to the effective photoelectric opacity 
actually occurs. As a function of H i column, a single power-law with slope —1.5 appears to provide 
a surprisingly good description over nearly 10 decades in iVni (Hu et al. 1995; Petitjean et al. 1993; 
Tytler 1987). In a recent analysis of the density distribution by redshift interval, Kim et al. (1997) 
have used high resolution, high S/N Keck spectra to show that the deficit of lines (relative to this 
fit) noted by Petitjean et al. above A^hi > 10^^'^ cm~^, while more pronounced at lower redshifts, 
tends to disappear (or moves to columns higher than 10-*^^ cm^^) at z > 3.5. At high redshift, it is 
then a good approximation to use for the distribution of absorbers along the line of sight: 



(5) 



- 5 - 



While it has been assumed in the past that the Lya forest clouds and LLS evolve at different 
rates (hence that there is an actual mismatch around 10^''cm~^ in the adopted function), recent 
determinations have derived a value of 7 = 1.55 it 0.39 for the LLS (Stengler-Larrea et al. 1995), 
consistent within the errors with 7 = 2.46ib0.37 obtained by Press &: Rybicki (1993), 7 = 1.89ib0.28 
by Bechtold (1994), and 7 = 2.78 ±0.71 by Kim et al. (1997) for the Lya forest line. For simphcity, 
we will assume here a single redshift exponent, 7 = 2, for the entire range in column densities. A 
normalization value of A'o = 4.0 x 10^ produces about 3 LLS per unit redshift at z = 3, as observed 
by Stengler-Larrea et al. (1995), and, at the same epoch, ~ 150 lines above A'^hi = lO^'^'^^ cm~^, in 
good agreement with the estimates of Kim et al. (1997). With this normalization and 7 = 2, the 
adopted distribution provides about the same H i photoelectic opacity as in the model discussed 
in Paper II, but underpredicts the attenuation of quasar continua due to Lya line blanketing (cf 
Paper I). 



2.3. Attenuation Length 

If we extrapolate the N^^'^ power-law in equation (|5|) to very small and large columns, the 
effective optical depth becomes an analytical function of redshift and wavelength, 
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where ao is the hydrogen photoionization cross-section at the Lyman edge ul, and we have not 
included the contribution of helium to the attenuation along the line of sight. [While the opacity 
due to He i is negligible in the case of a QSO-dominated background. He 11 absorption on the way 
must be included for Uo > 4z^i(l + Zo)/{l + z).] 

It is practical, for the present discussion, to define a redshift 2:1(2^1) such that the effective 
optical depth between Zq and zi is unity. When < i'l, photon emitted at zi with frequency 
J^i = 1^0(1 + •2i)/(l + Zo) will be redshifted below threshold before being significantly attenuated by 
intervening H I . From equation (|6|) it can be shown that the universe will be optically thin below 

^^^^^^=[(.lM)3/^-l]V3-^- 

This expression gives zi = 1.85 for vi = 1.2i/l and zi = 1.64 for vi = 2vl. It has a minimum for 
vi = 1.591/^, corresponding to what we shall term in the following the "breakthrough epoch" of the 
universe, Zbr = zi{l.bS)vL) = 1.56. For z < z^r, all radiation sources will be able to see each other 
in the hydrogen Lyman continuum. Hence, it is only at z < z^j. that the degree of ionization of the 
IGM will be determined by the balance between radiative recombinations and the total ionizing 
flux emitted (say) by all QSOs which appear after Zbr- 

By contrast, due to the rapid increase with lookback time of the number of absorbers, the 
mean free path of photons at 912 A becomes so small beyond a redshift of 2 that the radiation 
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is largely "local", as sources at higher redshifts are severely attenuated. Expanding equation (|6|) 
around z = 3, for example, one gets Teff (z^l) ~ 0.36(1 + z)'^Az = 1 for Az = 0.18. This corresponds 
to a proper distance or "absorption length" of only 

Al{uL)^33Mpc(^^^ (8) 

shorter than the mean free path between the rarer LLS because of the additional absorption from the 
numerous Lya forest clouds]^ The small absorption length is mostly due to systems with continuum 
optical depth around unity. In the local (or "source- function" ) solution to the equation of radiative 
transfer, this strong attenuation effect can be approximated by simply ignoring sources with z > zi 
and neglecting absorption for those sources with z < zi. Since only radiation sources within the 
volume defined by an absorption length contribute then to the background intensity shortward of 
the Lyman edge, cosmological effects such as source evolution and frequency shifts can be neglected, 
and it is possible to write AttJ^v) ~ e{y)/K,{u) = e{u)Al{u). In this approximation the number of 
ionizing photons per unit proper volume present at redshift z is given by 

riion = — [ dv^^^ = nion-^^, (9) 
c Ju]^ c 

where nion(0 = I!^ due{v,t) /hv^ and (A/) is an average over the incident photon spectrum.^ From 
the equation of ionization equilibrium, n\oa = ^H/^rec, an approximate relation can be derived 
between the absorption length and the volume-averaged gas recombination time i^ec if all absorbers 
are optically thin, highly ionized in both H and He, and contain most of the baryons of the universe, 

(AO^^cW, (10) 
where nn is the mean hydrogen density of the expanding IGM, nH(0) = 1.7 x 10~^ (175/150/0.08), 

i-ec = [(l + 2x)nHaBC]-i = 0.3Gyr(^^^ '(^)''c^3oS (H) 

as is the recombination coefficient to the excited states of hydrogen, x the helium to hydrogen 
cosmic abundance ratio, C = ('^hii)/'^hii ionized hydrogen clumping factor, ^ and a gas 

temperature of 10^ K has been assumed. Clumps which are dense and thick enough to be self- 
shielded from UV radiation will stay neutral and will not contribute to the recombination rate. 



^As filtering through a clumpy IGM significantly steepens the UV background spectrum (Miralda-Escude & 
Ostriker 1990; Madau 1992), the absorption length at the He 11 edge is smaller than at 1 ryd, AI^Avl) ~ 

1.5AZ(i/l)^J228/J912. 

^Since cx {v/vLf-^, one has (A?) = Al{uL)3{'i^'^-°'- - 1)/(1 - 4""' ), where 5 = a,/{1.5-as), e{u) oc 

and we have assumed a cutoff at 4 ryd because of He 11 absorption. A spectrum with = 1.8 yields (Al) = 2.2AI{vl)- 

^This may be somewhat lower than the total gas clumping factor if higher density regions are less ionized (Gnedin 
& Ostriker 1997). 
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An empirical determination of the dumpiness of the IGM at high redshifts is hampered by our 
poor knowledge of the ionizing background intensity and the typical size and geometry of the 
absorbers. Numerical N-body/hydrodynamics simulations of structure formation in the IGM within 
the framework of CDM dominated cosmologies (Cen et al. 1994; Zhang et al. 1995; Hernquist et 
al. 1996) have recently provided a definite picture for the origin of intervening absorption systems, 
one of an interconnected network of sheets and filaments, with virialized systems located at their 
points of intersection. In the simulations of Gnedin &; Ostriker (1997), for example, the clumping 
factor rises above unity when the collapsed fraction of baryons becomes non negligible, i.e. z < 20, 
and grows to (7 > 10 (40) at z ~ 8 (5) (because of finite resolution effects, numerical simulations 
will actually underestimate clumping): the recombination timescale is much shorter than that for 
a uniform IGM, and always shorter than the expansion time. 



2.4. Sources of Ionizing Radiation at High Redshifts 

24.1. The Quasar Cutoff 

The existence of a decline in the space density of bright quasars at redshifts beyond ~ 3 was 
first suggested by Osmer (1982), and has been since then the subject of a long-standing debate. 
In recent years, several optical surveys have consistently provided new evidence for a turnover in 
the QSO counts. The six-band multicolor survey of Warren, Hewett, & Osmer (1994, hereafter 
WHO) shows a drop by a factor of e~^-^^(^~^-^) beyond z 3.3. A decline by a factor of 2.7 per 
unit rcdshift in the range 2.75 < z < 4.75 is seen for Mb < —26 by the Schmidt, Schneider, &: 
Gunn (1995) grism survey. Kenncfick, Djorgovski, & de Carvalho (1995, hereafter KDC) found a 
decrease in the space density of QSOs with Mb < —27 by a factor of ~ 7 from z ~ 2 to z ~ 4.35, 
and no evidence for any luminosity dependent evolution. According to Hawkins & Veron (1996), 
the space density relation is rather flat in the range 1.5 < z < 3.0, and then gradually decline to 
z 4.5 for Mr < —25.5. The interpretation of the drop-off observed in optically selected samples 
is equivocal, however, because of the possible bias introduced by dust obscuration arising from 
intervening systems (Ostriker & Hcislcr 1984; Wright 1990; Fall & Pci 1993). Radio emission, on 
the other hand, is unaffected by dust, and Shaver et al. (1996) have recently shown that the space 
density of radio-loud quasars also decreases strongly for z > 3, demonstrating that the turnover is 
indeed real and that dust along the line of sight has a minimal effect on optically-selected QSOs. 

The comoving space density of bright quasars as a function of redshift is plotted in Figure 1. 
The data points are taken from Hartwick & Shade (1990) (a compilation of 15 major optical surveys 
with more than 1000 quasars covering the range 0.1 < 2; < 3.3), WHO (who have combined their 
own and four other surveys to form a sample which contains about 200 quasars at 2 < 2; < 4.5), 
Schmidt ct al. (1995) (90 quasars in the interval 2.75 ^ z < 4.75), and KDC (10 quasars at z > 4). 
All of them refer to QSOs with Mb < —26 except for KDC, whose survey limit is one magnitue 
brighter. The points have been normalized to the z = 2.5 space density of Mb < —26 quasars 
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{Mb < —27 in the case of KDC) as given by WHO. For comparison, similarly normalized space 
densities are shown for the Parkes fiat-spectrum radio-loud sample of Hook, Shaver, &; McMahon 
(1998). 

A successful analytical representations for the shape of the quasar blue luminosity function 
(LF) is the standard double power-law (Boyle, Shanks, &; Peterson 1988; Pei 1995): 

= wmMWim^' ^^^^ 

The data in Figure 1 can be well fit by a simple model in which the entire LF shifts along the 
luminosity axis as the position of the break evolves with redshift, 

L.(.)=L.(0)(l + .)---l^-^, (13) 

where a power-law spectral distribution for the typical quasar spectrum has been assumed, L{i') oc 
and the dependence on the spectral index as is explicitly shown. The comoving density of 
quasars brighter than an absolute magnitude limit as a function of redshift is given by 

N{z,MB<M,i^)= dMB(f>{MB,z). (14) 

In a cosmology with (/i5o,f/o) = (1)0.5), the LF fitting parameters values given by Pei (1995) are 
log (0*/Gpc-3) = 2.95, /3i = 1.64, (32 = 3.52, M*(0) = -22.35. The combination of exponentials 
we choose in this paper (to be compared with Pel's representation, a simple Gaussian form that 
appears to underestimate the number of quasars at z ^ 4.5) to describe the evolution of quasars 
over the whole range of redshifts has fitting parameters z^, = 1.9, ( = 2.58, and ^ = 3.16. The 
results of the model are plotted in Figure 1 for = 0.5. The density of QSOs has a relatively flat 
maximum at 1.8 ^ z < 2.8, and declines gradually at higher redshits. Note that, while at z < 2.5 
and z > 4 all estimates are in good agreement with each other, in the range 2.5 < z < 3.5 the 
Hartwick &; Shade and Schmidt et al. data points appear to be sistematically lower than the WHO 
point, a discrepancy that may reflects the uncertainties arising in the incompleteness corrections, 
different selection criteria, etc. 



2.4-2. Quasar Emissivity 

The QSO emission rate of hydrogen ionizing photons per unit comoving volume can then be 
written as 



(15) 



AfQ{z) = (1 + .)-- ^y + f)^Qi-B,0) _ / 
' ^ > e«^* ash ^ '\ 

where eqiyB^^) is the extrapolated 2 = emissivity in the blue band, 

/•QO 

eQ{i^B,0)= ^(L,0)L(iL « 9.6 X 10^^ ergs s-^Hz-^Mpc-^ (16) 

»' -^'min 
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and, as usual, we have assumed a cutoff in the spectrum at 4 ryd because of He ii absorption. The 
LF of Seyfert galaxies matches remarkably well that of optically selected QSOs at Mb = —23, and 
does not show clear evidence of leveling off down to ~ — 18 (Cheng et al. 1985). We therefore 
integrate the LF from Lmin ~ 0.018L,= (0), keeping in mind that, for the adopted shape of the LF, 
a value of Lmin 2 times smaller would only increase the quasar emissivity by about 5%. On the 
other hand, a steeper LF with /?i = 2.0 and same Lmin would boost eQ{i^B,0) up by a factor 1.75. 



Equation (15) assumes a single power-law from optical to EUV frequencies, a poor approx- 
imation to the "typical" quasar spectrum. In the following we shall adopt, except when stated 
otherwise, a more realistic form for the quasar spectral energy distribution: 

(2500 < A < 4400 A); 
1^-0-8 (1050 < A < 2500 A); (17) 
(A < 1050 A), 

where the different slopes have been continuously matched. This is based on the rest-frame optical 
and UV spectra of Francis et al. (1991) and Sargent, Steidel, & Boksenberg (1989), and is similar 
to the model used in Paper II except for a steeper index at the shortest wavelengths, a choice 
more consistent with the recent Hubble Space Telescope (HST) observations of the EUV spectra of 
radio-quiet QSOs at intermediate redshifts (Zheng et al. 1997; see also Laor et al. 1997). 

In Figure 2 we plot the value of A/g as a function of redshift according to what will become in 
the following our "standard QSO model". It is important to notice that the procedure adopted to 
derive this quantity implies a large correction for incompleteness at high-z. With a fit to the quasar 
LF which goes as (piL) oc L~^'^^ at the faint end, the contribution to the emissivity converges rather 
slowly, as LP'^^. At z = 4, for example, the blue magnitude at the break of the LF is « —25.4, 
comparable or slightly fainter than the limit of current high-z QSO surveys. A large fraction, 
about 90% at z = 4 and even higher at earlier epochs, of the ionizing emissivity in our model is 
therefore produced by quasars that have not been actually observed, and are assumed to be present 
based on an extrapolation from lower redshifts. While the value of Nq obtained by including the 
contribution from observed quasars only would be much smaller at high redshifts than shown in 
Figure 2, it is also fair to ask whether an excess of low-luminosity QSOs, relative to the best-fit 
LF, could actually boost the estimated ionizing emissivity at early epochs. The interest in models 
where the quasar LF significantly steepens with lookback time, and therefore predict many more 
QSOs at faint magnitudes than the extrapolation of Pel's (1995) fitting formulae, stems from recent 
claims of a strong linear correlation between bulge and observed black hole masses (Magorrian et al. 
1997), linked to the steep mass function of dark matter haloes predicted by hierarchical cosmogonies 
(Haehnelt & Rees 1993; Haiman & Loeb 1998; Haehnelt, Natarajan, & Rees 1998). We shall see 
in the next section how the space density of low-luminosity quasars at high-z may be constrained 
by the observed lack of red, unresolved faint objects in the Hubble Deep Field (HDF). 
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2.4-3. Quasar Candidates in the Hubble Deep Field 

Several surveys of unresolved objects in the HDF have appeared in the literature (Elson, 
Santiago, & Gilmore 1996; Flynn, Gould, & Bahcall 1996; Mendez et al. 1996). While the smah 
size of the field (about 5.3 square arcmin covered by the three WFC chips) does not allow one to 
set meaningful constraints to the number density of bright QSOs, its depth and resolution can be 
used to explore quasars counts to very faint magnitudes (Conti et al. 1998). In the following we 
shall denote with Bab, Vab, and Jab the AB magnitudes in the F450W, F606W, and F814W HDF 
bandpasses. 

Elson et al. find 59 unresolved objects down to a limiting magnitude of Iab ^ 27. 2.^ The bulk 
of them are faint and blue, {V — /)ab < 1.6 mag, with only the 9 brighter sources scattering out 
to {V — I)ab < 3.1. The number of point-like or nearly so objects identified by Mendez et al. 
and Flynn et al. is about 4 times smaller. Mendez et al. emphasize how, while the star-galaxy 
separation software is only reliable to Vab ~ 27.3, misclassification at fainter magnitudes turns 
more galaxies into stars than viceversa, as diffuse and faint features disappear first. The number 
of point-like sources detected in the range 27.3 < Vab < 29.8 should then be considered as an 
absolute upper limit to the number of stars (or quasars) observed in this magnitude interval. There 
are, however, no obvious candidate high-z QSOs in the published catalogues of unresolved objects 
in the HDF, as first remarked by Elson et al. In analogy with the case of star-forming galaxies, 
in fact, Lyman continuum absorption arising because of the ubiquitous presence of H i along the 
line of sight is so severe that faint quasars will drop out of the F450W image altogether by z ~ 4 
(Paper I; Haiman, Madau, & Loeb 1998). At z ^ 4.5, line blanketing from the Lya forest will also 
cause a significant apparent depression (AVab ~ 1 niag) in the F606W continuum image. High 
redshifts (z > 4) quasars are expected then to have very red B — V and reddish V — I colors. 
By contrast, all of the faint (/ab < 25.5) stellar sources in the Mendez et al. (1996) sample have 
{B — V)ab < 0.6 and {V — I)ab < 1-8. A recent detailed search for faint QSOs in the HDF by Conti 
et al. (1998) confirms the lack of a very red compact population: down to a 50% completeness limit 
of Vab = 29.6 {Iab = 28.6), no z > 4 quasar candidates are found. 

While a detailed comparison between deep imaging data and recent theoretical modeling of 
the evolution of the quasar LF has appeared elsewhere (Haiman et al. 1998), it is important in this 
context to take note of the interesting constraints that the lack of red, unresolved objects in the 
HDF may pose on reionization scenarios which rely on a large number density of low-luminosity 
quasars at early epochs. Figure 3 shows the predicted number-/AB magnitude relation of active 
galaxies in the redshift range 4 < z < 5.5 for the best-fit model discussed in the previous section 
and an illustrative variant, one that is characterized by a steeper LF with slope /?i = 2.0 and a 
comoving space density that remains constant above z = 2.5 instead of dropping. This alternative 
evolution model has been choosen in order to boost the actual emission rate of hydrogen-ionizing 



^Throughout this paper we shall denote with -Bab, Vab, and Iab the AB magnitudes in the F450W, F606W, and 
F814W HDF bandpasses. 
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photons from QSOs at 2; ~ 5 by about a factor of 5. It is clear from Figure 3 that this model predicts 
too many red faint stellar sources in the field of view. On the other hand, a large population of 
faint AGNs at high- 2; would be consistent with the data if, at these faint magnitude levels and high 
image resolution, the host galaxies of active nuclei could actually be resolved. 

We shall see in § 3 below how, taken at face value, our results on the production of photons 
above 1 ryd by QSOs may pose some problems to the idea that a quasar-dominated background 
is responsible for maintaining the intergalactic gas in a highly ionized state at early epochs. Here 
we will simply note that the ability of quasars to overcome recombinations depends sensitively on 
the dumpiness of the medium. The emission rate of ionizing photons in our standard QSO model, 
Nq PS 4x 10^° s"^ Mpc-3 at z ~ 5, is just a fraction of the total number of hydrogen recombinations 
in a clumpy, fast recombining universe, 72 1.5 x 10^^(r2fe/i5Q/0.08)^C3o s^^ Mpc~^ if T 10^ K. 

2.4-4- Star-forming Galaxies 

Galaxies with ongoing star- formation are another obvious source of Lyman continuum photons. 
The tremendous progress in our understanding of faint galaxy data made possible by the recent 
identification of star-forming galaxies at 2 < 2; < 4 in ground-based surveys (Steidel et al. 1996a) and 
in the HDF (Steidel et al. 1996b; Madau et al. 1996; Lowenthal et al. 1997) has provided new clues 
to the long-standing issue of whether galaxies at high redshifts can provide a significant contribution 
to the ionizing background flux. Since the rest-frame UV continuum at 1500 A (redshifted into 
the visible band for a source at z « 3) is dominated by the same short-lived, massive stars which 
are responsible for the emission of photons shortward of the Lyman edge, the needed conversion 
factor, about one ionizing photon every 10 photons at 1500 A, is fairly insensitive to the assumed 
IMF and is independent of the galaxy history for t ^ 10^-^ yr.^ 

A composite ultraviolet luminosity function of Lyman-break galaxies at z ~ 3 has been recently 
derived by Dickinson (1998). It is based on spectroscopically (about 375 objects) and photomet- 
rically selected galaxies from the ground-based and HDF samples, and spans about a factor 50 in 
luminosity from the faint to the bright end. Because of the uncertanties that still remain in the 
rescaling of the HDF data points to the ground-based data, the comoving luminosity density at 
1500 A is estimated to vary within the range 1.6 to 3.5 X 10^6 ergs s-iHz~iMpc~3_ ^he "best 
guess" Schechter fit gives log ((/)*/Gpc-3) = 6.1, a = 1.38, and MXb(1500) = -20.95 (Dickinson 
1998), the magnitude at the "break" corresponding to a star-formation rate slightly in excess of 
lOM0yr~^ (Salpeter IMF). Figure 4 shows the Lyman continuum luminosity function of galaxies 
at z 3 (at all ages > 0.1 Gyr one has L(1500)/L(912) 6 for a Salpeter mass function and 
constant star formation rate, Bruzual & Chariot 1998), compared to the distribution of QSO lu- 



''Shortward of the Lyman edge, however, the differences in the predicted ionizing radiation from model atmospheres 
of hot stars can actually be significant (see, e.g., Chariot 1996). 
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minosities at the same redshift. The comoving ionizing emissivity due to Lyman-break galaxies 
is 4.2 lb 1.5 X 10^^ ergss"^ Hz~^ Mpc~^, between 2 and 4 times higher than the estimated quasar 
contribution at z = 3. 

This number neglects any correction for dust extinction and intrinsic H I absorption. While 
it has been pointed out by many authors (e.g., Mcurcr ct al. 1997; Pcttini ct al. 1997; Dickinson 
1998) that the colors of Lyman-break galaxies are redder than expected in the case of dust-free star- 
forming objects, the prescription for a "correct" de-reddening is still unknown at present (note that 
redder spectra may also results from an aging population or an IMF which is rather poor in massive 
stars). A Salpeter IMF, Eb-y = 0.1 model with SMC-type dust in a foreground screen, for example, 
has been found to reproduce quite well the rest-frame ultraviolet colors of the HDF "UV dropouts" 
(Madau, Pozzetti, & Dickinson 1998). In this model the color excess £'912-1500 = 1.64£'b_v is 
rather small and can be safely neglected in correcting from observed rest-frame far-UV to the 
Lyman edge. However, for typical dust-to-gas ratios, it is the H i associated with dust that would 
greatly reduce the flux of Lyman continuum photons able to escape into the intergalactic space. 
The data point plotted in Figure 2 assumes a value of /esc = 0.5 for the unknown fraction of ionizing 
photons which escapes the galaxy H I layers into the intergalactic medium (e.g., Madau & Shull 
1996; Giallongo, Fontana, &: Madau 1997). The possible existence of a numerous population of 
galaxies below the detection threshold, i.e. having star formation rates < 0.5 Mq yr~^, with a space 
density well in excess of that predicted by extrapolating to faint magnitudes the a = 1.38 best-fit 
Schechter function, will be discussed in § 4. 

The LF of Lyman-break galaxies at 2; > 4 is very uncertain. An analysis of the B-band 
dropouts in the HDF - candidate star-forming objects at 3.5 < z < 4.5 - seems to imply a decrease 
in the comoving UV galaxy emissivity by about a factor of 2 in the interval 2.75 < z < 4 (Madau et 
al. 1996, 1998), but the error bars are still rather large. Adopting a L(1500) to L(912) conversion 
factor of 6, we estimate a comoving ionizing emissivity of 2 x 10^^/csc ergs Hz~^ Mpc~^ at 
2; « 4. One should note that, while a population of highly reddened galaxies at high redshifts 
would be missed by the dropout color technique (which isolates sources that have blue colors in 
the optical and a sharp drop in the rest- frame UV), it seems unlikely that very dusty objects with 
/esc ^ 1, such as the recently discovered population of sub-mm sources (e.g. Barger et al. 1998) 
would contribute in any significant manner to the ionizing metagalactic flux. 

3. Reionization of the Universe 

In the previous section we have reviewed the basic theory of cosmological radiative transfer 
in a clumpy universe, demonstrated the local character of the Lyman continuum flux at early 
epochs, and discussed the role of bright quasars and star-forming galaxies as candidate sources of 
photoionization at high redshifts. In this and the following sections we will try to shed some light 
on the time-dependent reionization process, and introduce the concept of a critical emission rate 
of hydrogen-ionizing photons per unit cosmological volume. 
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In inhomogeneous reionization scenarios, the history of the transition from a neutral IGM to 

one that is almost fully ionized can be statistically described by the evolution with redshift of the 
volume filling factor or porosity Q{z) of H II , He II , and He III regions. The radiation emitted 
by spatially clustered stellar-like and quasar-like sources - the number densities and luminosities 
of which may change rapidly as a function of redshift - coupled with absorption processes in a 
medium with a time-varying clumping factor, all determine the complex topology of neutral and 
ionized zones in the universe. When Q << 1 and the radiation sources are randomly distributed, 
the ionized regions are spatially isolated, every UV photon is absorbed somewhere in the IGM, 
and the ionization process cannot be described as due to a statistically homogeneous radiation field 
(Arons &: Wingert 1972; Meiksin & Madau 1993). While the final size of the expanding ionized 
bubbles is only limited by the individual source lifetime, the timescale for full recombination is 
rather short and relict ionization zones may recombine faster than they can accumulate (cf Arons 
&: Wingert 1972). As Q grows, the crossing of ionization fronts becomes more and more common, 
and the neutral phase shrinks in size until the reionization process is completed at the "overlap" 
epoch, when every point in space is exposed to Lyman continuum radiation and Q = 1. 



3.1. Time Evolution of Expanding Cosmological H ii Regions 

When an isolated point source of ionizing radiation turns on, the ionized volume initially grows 
in size at a rate fixed by the emission of UV photons, and an ionization front separating the H ii and 
H I regions propagates into the neutral gas. Most photons travel freely in the ionized bubble, and 
are absorbed in a transition layer. The evolution of an expanding H ii region is governed by the 
equation 

^^-3W, = ^-Ji, (18) 

at nu tree 

(Shapiro &: Giroux 1987), where Vi is the proper volume of the ionized zone, A^ion is the number 
of ionizing photons emitted by the central source per unit time, and all other symbols have their 
usual meaning. Across the I-front the degree of ionization changes sharply on a distance of the 
order of the mean free path of an ionizing photon. When tree ^ the growth of the H ii region 

is slowed down by recombinations in the highly inhomogeneous IGM, and its evolution can be 
decoupled from the expansion of the universe. Just like in the static case, the ionized bubble will 
fill its time-varying Stromgren sphere after a few recombination timescales, 

y7 = ^^^^(l-e-*/*-). (19) 

While the volume that is actually ionized depends on the luminosity of the central source, the time 
it takes to produce an ionization-bounded region is only a function of tree- One should point out 
that the use of a volume-averaged clumping factor in the recombination timescale is only justified 
when the size of the H ii region is large compared to the scale of the clumping, so that the effect 
of many clumps (filaments) within the ionized volume can be averaged over. This will be a good 
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approximation either at late epochs, when the IGM is less dense and H ii zones have had time to 
grow, or at earlier epochs if the ionized bubbles are produced by very luminous sources like quasars 
or the stars within halos collapsing from very high-o" peaks. From equation (|5|), the mean free path 
between absorbers having neutral columns > A^hi is 0.8 Mpc h^Q^ [(l+2;)/6)]~^'^ cm"^)'''^: 
it is only on scales greater than this value that the clumping can be averaged over. (Systems with 
-^Hi ^ 10^^ cm^^ do not dominate the mass, Miralda-Escude et al. 1996; if the overdensity is 
an increasing function of A'^hi, this implies, a fortiori, that Iow-A'^hi absorbers do not contribute 
significantly to the volume-averaged recombination rate.) On smaller scales underdense regions 
are ionized first, and only afterwards the UV photons start to gradually penetrate into the higher 
density gas. 



With these caveats in mind, equation ( |lq ) approximately holds for every isolated source of 
ionizing photons in the IGM. The filling factor of H ii regions in the universe, QhiI; is then equal 
at any given instant t to the integral over cosmic time of the rate of ionizing photons emitted per 
hydrogen atom and unit cosmological volume by all radiation sources present at earlier epochs, 
Jq ''T'ionit')dt' /n}i{t'), minus the rate of radiative recombinations, Jq Qmi{t')dt' /trecit')- Differenti- 
ating one gets 

^QhII _ ^ion _ Qmi ^2Q-) 
dt fill tree 

It is this simple differential equation - and its equivalent for expanding helium zones - that statis- 
tically describes the transition from a neutral universe to a fully ionized one, independently, for a 
given emissivity, of the complex and possibly short-lived emission histories of individual radiation 
sources, e.g., on whether their comoving space density is constant or actually varies with cosmic 
time. In the case of a time-independent clumping factor, equation (pO[) has formal solution 



Qmi{t) = I dt''^ew{-j- + J^] , (21) 




with tj-cc oc t'^. Assuming a constant comoving emissivity, this can be rewritten as 

(22) 



Qmiit) "^^""^ 



where Ei(x) is the exponential integral function. Contrary to the static case, the H ii regions will 
always percolate in an expanding universe with constant comoving emissivity and t — > oo. At high 
redshifts, and for an IGM with C ^ 1, one can expand around t to find 

Qmi{t) ~ -^^Uec- (23) 

The porosity of ionized bubbles is then approximately given by the number of ionizing photons 
emitted per hydrogen atom in one recombination time. In other words, because of hydrogen 
recombinations, only a fraction Uec/t (~ a few per cent at z = 5) of the photons emitted above 1 
ryd is actually used to ionize new IGM material. The universe is completely reionized when Q = 1, 
i.e. when 

n^oritrec = ^H- (24) 
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While this last expression has been derived assuming a constant comoving ionizing emissivity and 
a time-independent clumping factor, it is also valid in the case nion and C do not vary rapidly 
over a timescale tree- Figure 5 shows the H ii filling factor resulting from numerical integration of 
equation (pO|), as a function of redshift for the QSO photoionization model described in § 2.4, and 
a clumpy IGM with (constant) C = 1,10,20 and 30. The transition from a neutral to a ionized 
universe takes place too late in this model, as late as z = 3.4 for C = 30, and never before z = 4.5 
even in the limiting case of a uniform medium. One should stress here the illustrative character 
of the curves depicted in Figure 5; in the real universe the dumpiness of the ionized hydrogen 
component will increase with cosmic time as more baryons collapse, and may also be a function of 
the UV photon production rate (see discussion below). 

3.2. The Case of Centrally Concentrated Clumps 

The very simple model described in the previous section, one where individual H ii regions 
expand at a rate that depends on nion and on a mean clumping factor assumed independent of the 
UV production rate, may not provide a good description of the late stages of reionization, when the 
fraction of baryons collapsed in virialized structures is non negligible, the filling factor of ionized 
material is close to unity, and a significant fraction of all hydrogen recombinations (the rate of 
which is proportional to the square of the mass times the clumping factor of the recombining gas) 
may occur in the gas surrounding localized "islands" of neutral material with a small volume filling 
factor, e.g., the Lyman- limit absorbers. 

Consider, as an illustrative example, a fully ionized diffuse IGM after complete overlapping 
has occurred, but with denser, partially neutral spheres embedded in it, each characterized by 
an isothermal (r^^) radial dependence of the gas density. Irradiated from the outside by the 
background radiation field generated by the ionizing sources, the neutral core of each clump will 
be delimited by a "Stromgren surface" with radius proportional to J~^/^, which will shrink as the 
UV metagalactic flux builds up with time. The cross section of the neutral cores, which determines 
the mean free path A/ of ionizing photons, scales then as J~^/^ (hence A/ oc J^/^). Since J cx 
UV emissivity x A/, in this simple model J and AZ will scale at a given epoch as the cube and 
square of the UV production rate at that epoch, respectively. On timescales where the number 
density of clumps stays approximately constant, the clumping factor scales as the inverse cube of 
the Stromgren radius, increasing therefore as the cube of the photon production rate. Spheres with 
different density profiles show a change of behaviour when the radial power-law slope (gas density 
oc r", with a = —2 in the case of isothermal spheres) becomes flatter than a = —3/2. For these 
less peaked density distributions, the total recombination rate converges towards small radii: the 
Stromgren surface suddenly moves inward when the UV flux rises above some threshold, allowing 
the photon mean free path to increase rapidly. (In the limiting case a = —3/2 the clumping factor 
becomes independent of the emission rate of ionizing photons.) 

The equations put forward in the previous section provide a good description of the reionization 
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process only when the clumps are not centrally concentrated. In this case the recombination rate 
depends on the filling factor Q and saturates when most of the volume is ionized; J increases 
steeply, by a large factor, when Q becomes ~ 1 and the H ii regions overlap. In the presence of 
dense clumps with steep density gradients, however, recombinations occur mainly just outside the 
Stromgren surfaces (of course the high-density regions may be filaments rather than spheres, but 
that does not change the picture), and there is a different relationship between J and ?i.ion 

compared 

to that predicted in a model where the clumping does not depend on the ionizing emissivity. The 
impact of Lyman-limit absorbers with steep density profiles on the reionization histories depicted 
in Figure 5 will depend on their baryonic content as well as on their volume filling factor, both of 
which are only poorly known at the moment. 



3.3. Delayed He ii Reionization 

Because of its higher ionization potential, the most abundant (by a factor ~ 100) absorbing 
ion in the universe is not H I but He ii . The importance of intergalactic helium in the context of 
this work stems from the possibility of detecting the effect of "incomplete" He ii reionization in the 
spectra of z ~ 3 quasars as, depending on the dumpiness of the IGM (Madau & Meiksin 1994) 
and on the spectrum of UV radiation sources (Miralda-Escude & Rees 1993), the photoionization 
of singly ionized helium may be delayed until much later than for H i . 

Since H i and He i do not absorb a significant fraction of hv > 54.4 eV photons, the problem 
of He II reionization can be decoupled from that of other ionizations, and the equivalent of equation 
(|20| ) for He iii expanding regions becomes 

^QhcIII _ "-ion,4 _ QHelll ^25^) 
dt nHe irec,HeIII ' 

where nion,4 now takes into account only photons above 4 ryd, and trec.Heiii is ~ 6.5 times shorter 
than the hydrogen recombination timescale if ionized hydrogen and doubly ionized helium have 
similar clumping factors. It is interesting to note that, if the intrinsic spectrum of ionizing sources 
has slope L{ij) oc the first terms on the right-hand side of equations (p0|) and ( p5[ ) are 

actually equal, and a significant delay between the complete overlapping of H ii and He iii regions 
can only arise if recombinations are important. This effect is illustrated in Figure 5, where the 
expected evolution of the He iii filling factor is plotted for our standard QSO-photoionization model: 
He II reionization is never completed before z = 3 in models with C > 10. While a significant 
contribution to the UV background at 1 ryd from massive stars may push the hydrogen reionization 
epoch to higher redshifts, the ratio between the number of He ii and H i Lyman continuum photons 
emitted from star-forming galaxies is only about two percent (Leitherer & Heckman 1995), five 
times smaller than in typical QSO spectra. It is likely then that recent observations of intergalactic 
He II absorption at high redshift may provide a rather direct test of quasar evolution models. 

To date, a He ii absorption trough has been detected in the spectra of four distant QSOs. 
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The HST and Hopkins Ultraviolet Telescope observations do not resolve individual absorption lines, 
but rather provide an average optical depth below 304 A in the rest-frame: THeii = 1.00 ± 0.07 
towards HS1700+6416 [zq = 2.74, Davidsen, Kriss, & Zheng 1996), thcIi > 1.5 towards PKS 
1935-692 {zQ = 3.18, Jakobsen 1997), and thcII > 1.7 towards Q0302-003 {zq = 3.29, Jakobsen 
et al. 1994; Hogan, Anderson, & Rugers 1996). These, together with the spectrum of HE2347- 
4342 {zq = 2.89, Reimers et al. 1997), which reveals patchy absorption with low He ii opacity 
(rHeii ~ 0.5) "voids" alternating several Mpc sized regions with vanishing flux (thcII = ^-^2^)1 
suggest that helium absorption does not increase smoothly with lookback time, but rather in the 
abrupt manner expected in the final stages of inhomogeneous He 11 reionization by quasar sources. 
A detailed assessment of the implications of late He 11 reionization, which may also have received 
some support from the recent observation of a sharp drop in the ratio of Si iv to C iv absorption in 
the Lya forest sX z k. 2> (Songaila 1998), will be the subject of another paper. Here it is important 
to remark that, while delayed He 11 reionization in a clumpy IGM appears to be naturally linked 
to the observed decline in the space density of quasars beyond 2; ~ 3, the complete overlapping of 
He III regions occurs instead much earlier {z ^ 5) in models that predict many more faint QSOs 
at high redshifts (Haiman & Loeb 1998). 



4. Discussion 

We have seen in the previous sections that, in the approximation the clumping can be averaged 
over, only the photons emitted within one recombination timescales can actually be used to ionize 
new material. As tree <C t at high redshifts, it is possible to compute at any given epoch a critical 
value for the photon emission rate per unit cosmological comoving volume, A/iom independently of 
the (unknown) previous emission history of the universe: only rates above this value will provide 
enough UV photons to ionize the IGM by that epoch. One can then compare our determinations of 
A'ion to the estimated contribution from QSOs and star-forming galaxies. Equation (|^) can then 
be rewritten as 

The uncertainty on this critical rate is difficult to estimate, as it depends on the clumping factor of 
the IGM (scaled in the expression above to the value inferred at z = 5 from numerical simulations, 
Gnedin & Ostriker 1997) and the nucleosynthesis constrained baryon density. A quick exploration 
of the available parameter space indicates that the uncertainty on A/ion could easily be of order ±0.2 
in the log. The evolution of the critical rate as a function of redshift is plotted in Figure 2. While 
A/ion is comparable to the quasar contribution at z > 3, there is some indication of a significant 
deficit of Lyman continuum photons at z = 5. For bright, massive galaxies to produce enough 
UV radiation at 2; = 5, their space density would have to be comparable to the one observed at 
z ~ 3, with most ionizing photons being able to escape freely from the regions of star formation 
into the IGM. This scenario may be in conffict with direct observations of local starbursts below the 
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Lyman limit showing that at most a few percent of the stehar ionizing radiation produced by these 
luminous sources actually escapes into the IGM (Leitherer et al. 1995).^ If, on the other hand, faint 
QSOs with (say) Mab = —19 at rest-frame ultraviolet frequencies were to provide all the required 
ionizing flux, their comoving space density would be such (0.0015 Mpc~^) that about 50 of them 
would expected in the HDF down to Jab = 27.2. At z ^ 5, they would appear very red in y — I as 
the Lya forest is shifted into the visible. This simple model can be ruled out, however, as there is 
only a handful (7) of sources in the HDF with (V — I)ab > 1-5 mag down to this magnitude limit. 

It is interesting to convert the derived value of A/ion into a "minimum" star formation rate per 
unit (comoving) volume, (hereafter we assume Qbhlo = 0.08 and C = 30): 

-53.1 f-1 ^ n ni Q f-l 1^ + ^ 



P*{z) = AAion(^) X 10-5^- Vc^c « 0.013/,-,^ (^-^ j Mq yr-1 Mpc'^. (27) 

The conversion factor assumes a Salpeter IMF with solar metallicity, and has been computed using 
Bruzual & Chariot (1998) population synthesis code. It can be understood by noting that, for 
each 1 Mq of stars formed, 8% goes into massive stars with M > 2OM0 that dominate the Lyman 
continuum luminosity of a stellar population. At the end of the C-burning phase, roughly half of 
the initial mass is converted into helium and carbon, with a mass fraction released as radiation of 
0.007. About 25% of the energy radiated away goes into ionizing photons of mean energy 20 eV. 
For each 1 Mq of stars formed every year, we then expect 

0.08 X 0.5 X 0.007 x 0.25 x Mqc^ 1 „ , 1 

— ~ 10" phots- (28) 



to be emitted shortward of 1 ryd. Note that the star formation density given in equation 
is comparable with the value directly "observed" (i.e., uncorrected for dust reddening) at z ~ 3 
(Madau et al. 1998; Dickinson 1998). If the volume emissivity has spectrum e(z^) oc i^""", the 
"minimum" background intensity at the Lyman limit is, from equations (^), (^, and ([2^ ) 

Jguiz) = ^Mon{z)has{l + zfAl{uL) - W-^^ as (^^^^ ^ ergscm'^ g-^ Hz-^ sr-\ (29) 

The same massive stars that dominate the Lyman continuum flux also manufacture and return 
most of the metals to the ISM. In the approximation of instantaneous recycling, the rate of ejection 
of newly sinthesized heavy elements which is required to keep the universe ionized at redshift z is, 
from equation (p7|), 

pz{z) = y{l - R)p.{z) > 3.5 x 10"^ (^^) (^i±f /^-i Meyr-^ Mpc'^, (30) 

where y is the net, IMF-averaged "yield" of returned metals, Zq = 0.02, and R ~ 0.3 is the mass 
fraction of a generation of stars that is returned to the interstellar medium. At z = 5, and over 



6 



Note that, at 2: = 3, Lyman-break galaxies would radiate more ionizing photons than QSOs for /esc ^ 30%. 
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a timescale of At = 0.5 Gyr corresponding to a formation redshift Zf = 10, such a rate would 
generate a mean metallicity per baryon in the universe of 



comparable with the level of enrichment observed in the Lya forest at z 3 (Songaila 1997): more 
than 2% of the present-day stars would need to have formed by 2; ~ 5. It has been suggested (e.g., 
Miralda-Escude k. Rees 1998) that a large number of low-mass galactic halos, expected to form 
at early times in hierarchical clustering models, might be responsible for photoionizing the IGM 
at these epochs. According to the spherically-symmetric simulations of Thoul & Weinberg (1996), 
photoionization heating by the UV background flux that builds up after the overlapping epoch 
completely suppresses the cooling and collapse of gas inside the shallow potential wells of halos 
with circular velocities ^ 35 kms~^. Halos with circular speed Vdrc = 50 kms~^, corresponding 
in top-hat spherical collapse to a virial temperature T^i^- = ^.h^mpv'^^^^/k ~ 10^ K and halo mass 
M = Q.lvl^JGH « 4 X 109[(1 + z)/Q]-^l^hv2 Mq, appear instead largely immune to this external 
feedback (but see Navarro & Steinmetz 1997). In these systems rapid cooling by atomic hydrogen 
can then take place and a significant fraction, /fift, of their total mass may be converted into stars 
over a timescale At comparable to the Hubble time (with / close to 1 if the efficiency of forming 
stars is high). If high-z dwarfs with star formation rates fQhM/At ~ 0. 3/0.5 Atj^g yr~^ were 
actually responsible for keeping the universe ionized Sit z ^ 5, their comoving space density would 
have to be 



two hundred times larger than the space density of present-day galaxies brighter than L*(4400), and 
about five hundred times larger than that of Lyman-break objects at 2; fa 3 with M < M^3(1500), 
i.e. with star formation rates in excess of 10 Moyr^^. Only a rather steep luminosity function, 
with Schechter slope a ~ 2, would be consistent with such a large space density of faint dwarfs and, 
at the same time, with the paucity of brighter B- and F-band dropouts observed in the HDF. The 
number density on the sky would be ^ 0.2 arcsec~^, corresponding to more than three thousands 
sources in the HDF. With a typical apparent magnitude at 2; = 5 of Iab ^ 29.5 mag (assuming 
/ = 0.5), these might be too faint to be detected by the HST, but within the range of the proposed 
Next Generation Space Telescope (Stockman et al. 1998). 

These estimates are in reasonable agreement (although towards the low side) with the calcula- 
tions of Miralda-Escude & Rees (1998), who used the observed metallicity of the Lya forest as the 
starting point of their investigation. A higher density of sources - which would therefore have to 
originate from lower amplitude peaks - would be required if the typical efficiency of star formation 
and/or the escape fraction of ionizing photons were low, (/, /esc) ^ 1- In this case the dwarfs could 
still be detectable if a small fraction of the gas turned into stars in very short bursts (there would 
then be an extra parameter associated with their duty cycle, in addition to /esc and /). A reduction 
of the star formation rate in halos with low circular velocities (necessary in hierarchical cosmogonies 




(31) 




(32) 
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to prevent too many baryons from turning into stars as soon as the first levels of the hierarchy 
collapse, White & Prenk 1991) may result from the heating and possible expulsion of the gas due 
to repeated supernova (SN) explosions after only a very small number of stars have formed. Recent 
numerical simulations (Mac Low & Ferrara 1998) show, however, that metal-enriched material from 
SN ejecta is accelerated to velocities larger than the escape speed from such systems far more easily 
than the ambient gas. If the same population of dwarf galaxies that may keep the universe ionized 
at z ^ 5 were also responsible for polluting the IGM with heavy elements, it is interesting to ask 
whether these atoms could diffuse uniformly enough to account for the observations of weak but 
measurable C IV absorption lines associated with the Lyct forest clouds (Lu et al. 1998; Songaila 
1997). Our fiducial estimate of 0.1 sources per comoving Mpc'^ implies that, in order to pollute the 
entire IGM, the metal-enriched material would have to be expelled to typical distances of order 1 
Mpc, which would require an ejection speed of about 200[(1 + 2;)/6]^/^ kms~^. In fact, the heavy 
elements may be restricted to filaments within only (say) 20 proper kpc (100-150 comoving kpc) 
of the halos, and in this case it would be enough to accelerate the outflowing metals to velocities 
comparable to the escape velocity, rather than to the higher speeds associated with SN ejecta. The 
required diff'usion of heavy elements would be a more serious constraint if the relevant galaxies were 
rarer and more luminous, as would happen if they originated from 3-cr peaks and star formation 
was postulated to occur with higher efficiency than in more typical peaks. 

We have benefited from useful discussions with A. Conti, M. Haehnelt, J. Miralda-Escude, and 
especially D. Sciama. Support for this work was provided by NASA through ATP grant NAG5-4236 
(PM), and by the Royal Society (MJR). 
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Fig. 1. — Comoving space density of bright QSOs as a function of redshift for (/150, Q'O) cts) = 
(1.0, 0.5, 0.5). The data points with error bars are taken from Hartwick &; Shade (1990) (filled dots), 
WHO (filled squares), Schmidt et al. (1995) (crosses), and KDC (filled pentagon). The points have 
been normalized to the z = 2.5 space density of quasars with Mb < M^^n = —26 {Mb < Mny^ = —27 
in the case of KDC) as given by WHO. The curves {solid line for M\[^ = —26, and dotted line for 
Miim = -27) have been computed from eqs. (|l|), (|l|), and (0), with log {(j)^/Gpc-^) = 2.95, 
Pi = 1.64, (32 = 3.52, M*(0) = -22.35, = 1.9, C = 2.58, and ^ = 3.16, and have been normalized 
to WHO as well. The empty triangles show the space density (normalized to the peak) of the 
Parkes flat-spectrum radio-loud quasars with P > 7.2 x 10^^ W Hz~^ sr~^ (Hook et al. 1998). 
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Fig. 2. — Comoving emission rate of hydrogen Lyman continuum photons {solid line) from QSOs, 
compared with the minimum rate {dashed line) which is needed to fuhy ionize a fast recombining 
(with clumping factor C = 30) Einstein-de Sitter universe with h^Q = 1 and 0,h = 0.08. Models 
based on photoionization by quasar sources appear to fall short at z = 5. See text for details on the 
assumed quasar luminosity function and spectral energy distribution. The data point with error 
bars show the estimated contribution of star- forming galaxies at z « 3. The fraction of Lyman 
continuum photons which escapes the galaxy H i layers into the intergalactic medium is taken to 
be /esc = 0.5. 
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Fig. 3. — Cumulative number-magnitude relation of quasars in the redshift range 4.0 < z < 5.5. 
The solid line shows the prediction for our standard model, one in which the faint end of the QSO 
luminosity function has slope (3i = 1.64 and tracks the turnover observed in the space density of 
bright quasars at z > 3. The dashed line shows the increased number of sources expected in the 
case of a steeper, /5i = 2.0, luminosity function and a comoving space density which stays constant 
above z = 2.5. The latter evolution scenario provides, within the errors, enough UV photons to keep 
the universe ionized at 2; ~ 5, but might be inconsistent with the lack of red, faint stellar objects 
observed in the Hubble Deep Field. Note that, in the assumed cosmology, Jab = 28 corresponds 
to an absolute mag nitude of Mab(1360A)= -18.2 at z = 5, at the very faint end of the AGN 
luminosity function. 
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Fig. 4. — The 912 A luminosity function of galaxies at z ~ 3 {solid line), compared to the distribu- 
tion of QSO luminosities at the same redshift {dashed line). The latter has been derived assuming 
a spectral slope of = 0.5. The former assumes a Salpeter IMF with constant constant star for- 
mation rate (age=l Gyr): Mab(912A)= —19 corresponds to a rate of ISMQyr^"*^. The solid and 
dashed lines represent functional fits to the data points, and the dotted lines their extrapolation. 
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Fig. 5. — The evolution of the H ii (solid lines) and He iii [dashed lines) filling factors as a function 
of redshift for an inhomogeneous universe where photoionization is dominated by QSOs, and clumps 
with steep density gradients make only a negligible contribution to the total recombination rate. 
The IGM density is taken to be ^bh^Q = 0.08, and, from right to left, the four curves assume a 
time-independent clumping factor of C = 1, 10, 20, and 30 (C = 1 in the case of a uniform IGM). 
The QSO intrinsic spectrum varies as v^^'^ shortward of the hydrogen Lyman edge. Note how the 
ionization of H I (He ii ) is never completed before z = 4.0 (z = 3.0) in models with C > 10. 



